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We calculate the one-body temperature Green's (Matsubara) function of the unitary Fermi gas 
via Quantum Monte Carlo, and extract the spectral weight function A(p^(jj) using the methods of 
maximum entropy and singular value decomposition. From A(p,a;) we determine the quasiparticle 
spectrum, which can be accurately parametrized by three functions of temperature: an effective 
mass m*, a mean-field potential [/, and a gap A. Below the critical temperature Tc = O.ISsf 
the results for m*, U and A can be accurately reproduced using an independent quasiparticle 
model. We find evidence of a pseudogap in the fermionic excitation spectrum for temperatures up 
to T* ^ 0.20£F > Tc. 
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Over the last few years, the field of cold atoms has 
drawn unprecedented attention, as documented in two 
recent review articles [1]. Within this field, the uni- 
tary Fermi gas, defined as the limit of vanishing inter- 
action range and infinite scattering length, continues to 
be a fascinating area of research for a number of reasons. 
First, the properties of the unitary regime are universal, 
making this problem relevant to a wide range of fields 
including string theories, the quark-gluon plasma, neu- 
tron stars, nuclei, and to a certain extent to high Tc- 
superconductors. Secondly, experimentalists can control 
the strength of the interaction by means of Feshbach res- 
onances, which allows for the systematic exploration of 
weakly as well as strongly coupled regimes. Finally, these 
systems exhibit a rich variety of phenomena and proper- 
ties (many of which await verification), creating an ideal 
playground for a large set of many-body techniques, pos- 
sibly the largest ever applied to a single problem. 

Properties established so far include: energy as a func- 
tion of temperature, entropy, frequencies of collective 
modes, speed of sound, critical temperature for the onset 
of superfluidity, and moment of inertia. In the case of 
polarized Fermi systems, the critical spin polarization at 
which superfluidity disappears has also been determined. 
All of these properties have been established quantita- 
tively with a reasonable degree of certainty and accuracy, 
both experimentally and theoretically. 

In spite of great efforts on the part of both theorists 
and experimentalists, some fundamental properties of 
these systems remain unknown. Among the most press- 
ing questions is the magnitude of the pairing gap and its 
evolution with temperature. Theoretical progress in this 
direction has been nearly at a standstill, except for the 
theoretical determination of the pairing gap at T = in 
Refs.jlJS, 0, m, the recent analysis of experimental data 
of Ref . and a recent experiment [t^ . 

While theoretical models abound, predictions are 



mostly qualitative, and their validity and accuracy are 
difficult to assess due to the absence of a small param- 
eter for a Fermi gas at unitarity. This work presents 
the first ah initio evaluation of the one-body temperature 
propagator of the unitary Fermi gas, free of uncontrolled 
approximations, which allows for the extraction of the 
temperature dependence of the pairing (pseudo)gap Q. 

We begin by defining the one-body temperature 
Green's (Matsubara) function 

g{p,T) = lTr{exp[-(/3 - t){H - iJiN)]i,^{p) x 

eM-T{H - ^JiN)i,{p)]], (1) 

where /3 = 1/T is the inverse temperature and r > 0. 
The trace Tr is performed over the Fock space, and 
Z = Tv{exY>[— (3{H — ^N)]} . The spectral weight function 
A{p, oj) can be extracted from the temperature Green's 
function using the relation: 
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By definition, A(p, uo) fulfills the following constraints: 

-^A{p,u;) = l. (3) 

Since our study focuses on the spin-symmetric system, 
and the Hamiltonian is spin-symmetric as well, G{p^ t) is 
diagonal in the spin variables and these are suppressed 
in all formulas. The numerical evaluation of the one- 
body temperature propagator ([1]) is performed as de- 
scribed in Refs.jlO, El, by using a Trotter expansion of 
exp[— r(i^ — /iA^)], followed by a Hubbard- Stratonovich 
transformation of the interaction and an evaluation of 
the emerging path-integral via Metropolis importance 
sampling. The number of imaginary time steps required 
to obtain an accuracy smaller than the statistical error 



2 



varies with temperature. At low temperatures the num- 
ber of time steps is O(IO^) [TH- All calculations pre- 
sented here have been performed with an average total 
particle number of 50-55 on an 8^ lattice with periodic 
boundary conditions [10]. We have generated between 
6000 and 10000 uncorrelated samples at each tempera- 
ture and the statistical errors are typically below 1%. 
The systematic errors, some due to finite lattice effects, 
others due to finite range effects, are estimated at about 
10-15%. Our T = extrapolation results [uj for the 
energy per particle are systematically lower than previ- 
ous fixed node Monte Carlo results which are variational 
[I, 0, [lH . We have not used the fixed- node approxima- 
tion and the value for ^ = 5E/3NeF ~ 0.40 that we 
extract at unitarity is in agreement with the auxiliary 
field Monte Carlo results of Ref. [iS^. 

The numerical determination of co) via inversion 
of Eq. ([2|) is an ill-posed problem that requires spe- 
cial methods. We have used two, based on completely 
different approaches. The first approach is the max- 
imum entropy method [14], which is based on Bayes' 
theorem. Quantum Monte Carlo (QMC) calculations 
provide us with a discrete set of values r^), where 
i = l,2,...,A/'r = 50. We treat them as normally dis- 
tributed random numbers around the true values Q{p^ Ti). 
The Bayesian strategy consists in maximizing the pos- 
terior probability P{A\G) ex P(G\A)P{A) of finding the 
right uo) under the condition that ^(p, Ti) are known. 
Here, P{G\A) ex exp(— ^x^) is the likelihood function^ 

where = Y.ti [Q{P.n) - g{p,n)^ / cj\ The quan- 
tity Qip^Ti) is determined by the spectral weight func- 
tion in the discretized form of Eq. ([2]) at frequencies 
uok' The prior probability describing our igno- 

rance about the spectral weight function, is defined as 
P{A) oc exp(a6'(A1)), where a > and S{M) is the rel- 
ative information entropy with respect to the assumed 
model Ai: 



The operator /C possesses a singular system defined as: 



S{M) = "^Acj 
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Hence the maximization of P{A\G) leads in practice to 
the minimization of the quantity ~ aS{M.) with re- 
spect to A. Note that the parameter a governs the rela- 
tive importance of the two terms. The entropy term pre- 
vents excessive inclusion of unjustified structure into the 
shape of the spectral weight function. The constraints 
([3]) are enforced by means of Lagrange multipliers. 

The second approach is based on the singular value 
decomposition of integral kernel JC of Eq. ([2|), which can 
be rewritten in operator form as 



(6) 



where /C* denotes the adjoint of /C, the are the singular 
values and the Ui^ Vi are right-singular functions and left- 
singular vectors respectively. The singular system forms 
a suitable basis for the expansion of the spectral weight 
function [is'], which we can then write as 

A{p,u;)=J2b^{p)ui{u;), = l(g(p) . (7) 
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where [ • ] is a scalar product and r is the rank of the op- 
erator /C/C*. Since Q{p^ Ti) is affected by the QMC errors 
(Ji, the coefficients hi carry some uncertainty Abi. Each 
set of expansion coefficients bi G {bi — Abi^ bi + Abi) re- 
produces Q{p^ Ti) within its error bars. We use this flexi- 
bility of choosing the expansion coeflicients to produce a 
solution satisfying constraints (|3]) [l6|. 

The advantages and disadvantages of both methods 
will be discussed elsewhere [17]. Here we note only that 
since they are based on completely different approaches 
their agreement serves as a robust test for the determi- 
nation of the spectral weight function. A sample of cal- 
culated spectral weight functions at unitarity are shown 
in Fig. [H In order to characterize the quasiparticle exci- 
tation spectrum we have associated with the maximum 
of A{p^uo) the quasiparticle energy E{p)\ 



E{p) = ±^ 
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Q{P,Ti) = {K,A){p,Ti). 



(5) 



where m* is an effective mass, the potential U and the 
"pairing" gap A depend on temperature, and fi is an in- 
put parameter. In Fig. [2] we compare the spectrum of 
elementary fermionic excitations evaluated in Ref. 
with the one extracted by us from our lowest tempera- 
ture spectral weight function. Such comparisons are le- 
gitimate because the temperature dependence of various 
quantities at T < Tc is relatively weak, see Ref. [lO] and 
the results below. The agreement between the T = 
Monte Carlo results and the low T-limit of our finite- 
T unrestricted QMC data validates those results. One 
should keep in mind that the effective range corrections 
are noticeable, since ro = 4/7rA ^ 0.4 (in lattice units), 
where A is the cutoff in momentum [10]. 

The extracted value of the single-particle potential U 
(see Fig. [3j) shows essentially no temperature dependence 
in the range investigated by us. (Simulations at higher 
temperatures are not reliable with lattice sizes we con- 
sidered here, see also Ref. [iq|-) The values of U and a 
are very close to the values determined in Ref. [18] at 
T = and show almost no temperature dependence. 

A surprising feature of our results can be seen if we as- 
sume that the system is composed of independent quasi- 
particles with BCS-like dispersion relation, in which case 
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FIG. 1: (Color online) Spectral weight function A{p,oo) for 
three temperatures: T = O.ISsf ~ Tc (upper panel), T = 
0.18s F (middle panel) and T = 0.20s f (lower panel). The 
presence of a gap in clearly seen in the upper two panels. 



the susceptibility can be easily evaluated: 



x{p) 



/3 



drQ{p, r) 



E{p) e/5^(p) + 1 • 



(9) 



From the calculated one-body propagator, using Eqs. 
()l|9p . one can extract the spectrum of the elementary 
fermionic excitations at finite temperature, which turns 



2.5 
2 

m 1 



^ This work: T=0.18|_ 
E(p) fit 

• Carlson and Reddy 




e(p)/ep 



FIG. 2: (Color online) Quasiparticle energies E{p) (squares) 
extracted from the spectral weight function A(p,uj) at T = 
0.1s F- The line corresponds to the fit to Eq. JH]). The circles 
are the results of Carlson and Reddy 01 • 
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FIG. 3: (Color online) The single-particle parameters ex- 
tracted from the spectral weight function at unitarity. The 
dashed, solid and dotted lines represent the quantities: m/m*, 
A/sf: U/sf, respectively, extracted using the assumption of 
independent quasiparticle model. 



out somewhat unexpectedly to be accurately parameter- 
ized by Eq. ([8]). Extracted parameters are plotted in 
Fig. [3] with lines. While the agreement between the 
mean-field potential U and effective mass m* obtained 
using the two procedures is almost perfect at all temper- 
atures, the pairing gap is reproduced satisfactorily only 
up to T < Tc. 

Another notable feature of our results is that both 
methods admit the "gapped" spectral function above the 
critical temperature, routinely referred to as the pseudo- 
gap [8] . Various aspects and the physics of a pseudogap in 
a Fermi gas in the unitary regime have been advocated 
and discussed for a number of years by several groups 
[Tot . It is however notable that the pseudogap has not 
entered the mainstream of research in this field and the 
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physics of the pseudogap is barely covered in the recent 
reviews [ij , which reflects a rather widespread opinion in 
the cold atom community that the gap should vanish at 
Tc- There have been several experimental attempts to ex- 
tract the pairing gap in ultracold dilute Fermi gases [2q| 
and a theoretical explanation of these spectra was given 
in Refs. [21|]. However, it was later shown in Refs. [22[ 
that these initial interpretations of the rf-spectra as re- 
vealing the pairing gap were in error, as strong final 
state interaction effects had been neglected. In partic- 
ular the difficulties of determining the pairing gap using 
rf-spectroscopy were discussed in Ref. [23|. In Ref. [2J] 
the use of Bragg spectroscopy is advocated instead, in 
order to measure the onset of superfluidity as well as the 
appearance of a pseudogap. A recent theoretical calcu- 
lation of the spectral function, based on summing all the 
particle-particle ring diagrams [25], does not reveal signs 
of a pseudogap. On the other hand in Refs. 2^, 2^ there 
are experimental indications that a pseudogap exists in 
a unitary gas. 

Our calculations show that the spectral function re- 
veals the presence of a gap in the spectrum up to about 
T* ^ 0.20£f. This result is reproduced by both the max- 
imum entropy and singular value decomposition meth- 
ods, which indicate that the spectral function possesses 
a two peak structure around the Fermi level at temper- 
atures above Tc. In the maximum entropy method this 
result is stable with respect to variation of the assumed 
model M. and for a^a < 0.3, see Eq. (|4]). In range of 
temperatures 0.18 . . . 0.25£f the singular value decom- 
position method can reproduce both two- and one-peak 
structures depending on the details of the method, due 
to finite size of the statistical errors in the imaginary 
time Green's function. In our case this resolution is at 
the level of A 0.2 . . . O.S^f, which means that the gap 
above Tc is barely visible by the singular value decom- 
position method. We note that T* is the temperature 
at which, not surprisingly in hindsight, the caloric curve 
E{T) has a shoulder [lO|, which in Ref. 
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